Quantum Arnol'd Diffusion in a Simple Nonlinear System 
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We study the fingerprint of the Arnol'd diffusion in a quantum system of two coupled nonlinear 
oscillators with a two-frequency external force. In the classical description, this peculiar diffusion 
is due to the onset of a weak chaos in a narrow stochastic layer near the separatrix of the coupling 
resonance. We have found that global dependence of the quantum diffusion coefficient on model 
parameters mimics, to some extent, the classical data. However, the quantum diffusion happens 
to be slower that the classical one. Another result is the dynamical localization that leads to a 
saturation of the diffusion after some characteristic time. We show that this effect has the same 
nature as for the studied earlier dynamical localization in the presence of global chaos. The quantum 
Arnol'd diffusion represents a new type of quantum dynamics and can be observed, for example, in 
2D semiconductor structures (quantum billiards) perturbed by time-periodic external fields. 



PACS numbers: 05.45-a, 03.65-w 



I. INTRODUCTION 



As is well known, dynamical chaos in Hamiltonian sys- 
tems is related to the destruction of separatrices of non- 
linear resonances. For a strong interaction between the 
resonances, the latter can overlap in the phase space, 
thus leading to a global chaos for which chaotic region is 
spanned over the whole phase space of a system, although 
large isolated islands of stability may persist. For a weak 
interaction, chaotic motion occurs only in the vicinity 
of separatrices of the resonances, in accordance with the 
Kolmogorov-Arnol'd-Moser (KAM) theory (see, for ex- 
ample, 0). 

In the case of two degrees of freedom (N — 2), a pas- 
sage of the trajectory from one stochastic region in phase 
space to another is blocked by KAM surfaces. The sit- 
uation changes drastically in many-dimensional systems 
(N > 2), where the KAM surfaces no longer separate 
one stochastic region from another, and chaotic layers 
of destroyed separatrices form a stochastic web that can 
cover the whole phase space. Thus, if trajectory starts 
inside the stochastic web, it can diffuse throughout the 
phase space. This weak diffusion along stochastic webs 
was predicted by Arnol'd in 1964 and since that time 
it is known as a very peculiar phenomenon, universal for 
nonlinear Hamiltonian systems with N > 2. 

There are many physical systems which behavior ap- 
pears to be strongly affected by the Arnol'd diffusion. As 
examples, one should mention three-boby gravitational 
systems || and galaxies dynamics Q. It is argued that 
the Arnol'd diffusion may have strong impact for the be- 
haviour of our solar system, it is also responsible for a 
loss of electrons in magnetic traps (see discussion and 
references in Q). From the practical point of view, the 
Arnol'd diffusion may be dangerous for long-time stabil- 
ity of motion of charge particles in high energy storage 
rings §. 

Recently, the Arnol'd diffusion was explored in the 



classical description of Rydberg atoms placed in crossed 
static electric and magnetic fields j7|. The semiclassical 
approach has been used for the stochastic pump model 
S , where the effect of quantization of the Arnol'd diffu- 
sion in a system of two pairs of weakly coupled oscillators 
has been investigated. 

For the first time, the Arnol'd diffusion was observed 
in numerical study of a 4-dimensional nonlinear map || . 
The more physical model of two coupled nonlinear oscil- 
lators with time-dependent perturbation was considered, 
both analytically and numerically, in (n]jll[] (see also re- 
view [El and the books [jlj). Numerical experiments with 
this model have confirmed analytical estimates obtained 
for the diffusion coefficient in dependence on model pa- 
rameters (for recent studies on this subject see p2[). 
Note that direct numerical study of the Arnol'd diffusion 
is quite difficult since its rate is exponentially small, and 
it occurs only for initial conditions inside vary narrow 
stochastic layers. 

One should distinguish stochastic web for the Arnol'd 
diffusion from that found for systems which are linear in 
the absence of a perturbation, see for example |fb3|| . In 
the latter case the stochastic web arises around a single 
nonlinear resonance that has infinite number of cells in 
the phase space (see Jl^j and references therein). 

So far, the studies of the Arnol'd diffusion are re- 
stricted by the classical (or semiclassical) approaches. 
On the other hand, it is important to understand the 
influence of quantum effects. The problem is not triv- 
ial since strong quantum effects can completely suppress 
weak diffusion along narrow stochastic layers, even in a 
deep semiclassical region The purpose of this work 
is to study the fingerprint of the Arnol'd diffusion in a 
quantum model, by making use a direct numerical sim- 
ulation of some nonlinear system, both in classical and 
quantum description. Preliminary data are reported in 
Ref. 0. 

The paper is organized as follows. In Sec. II the basic 
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model is introduced and discussed. We describe shortly 
the mechanism of the classical Arnol'd diffusion in clas- 
sical system and give the expression for the diffusion co- 
efficient. Sec. Ill is devoted to the study of the quan- 
tum model. First, in Sec. III(A) we study the structure 
of eigenstates of the stationary model (without an ex- 
ternal field). Second, in Sec. HI(B) we show how to 
construct the evolution operator that allows us to inves- 
tigate the dynamics of the system. Here we also discuss 
global properties of quasienergy states. Next step of our 
consideration (Sec. 111(C)) is the study of the evolution 
of the system for different initial conditions and model 
parameters. We show that for initial states correspond- 
ing to stochastic layer near the separatrix of the coupling 
resonance, the motion has a diffusion-like character. We 
calculate quantum diffusion coefficient and compare it 
with the classical one. Quantum effects of the dynamical 
localization and suppression of the classical diffusion are 
discussed in Sec. III(D). In Sec. IV we give our conclu- 
sions by summarizing main results, and shortly discuss 
possible systems in which the quantum Arnol'd diffusion 
may be observed. 
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where a is the amplitude of oscillations and K{l/\2) 
stands for the complete elliptic integral of the first kind. 
One can see that the amplitude a m of higher harmonics 
sharply decreases with an increase of m. Therefore, in 
action-angle variables I x , Q x one can approximately rep- 
resent the position of the x— oscillator by the expression 
x rs a(I x ) cos Q x . Then the system under consideration 
is described by the Hamiltonian 



H = A{lt 



4/3 



h ) (4) 
-fj,a[I x )a(Iy) cos Q x cosQ y - a(I x ) cos6 x /(t), 



with A = 



Near the coupling resonance uj x — uj y the resonance 
phase Q x — Q y and amplitudes a x , a y oscillate. Thus, it 
is convenient to introduce slow {6\ = <d x — & y ) and fast 
(9 2 = <d x + <d y ) phases by making use of the canonical 
transformation with the generating function, 



F = (Q x — & v )h + (Q x + Q y )h 



(5) 



II. CLASSICAL MODEL 

In this Section we discuss main results on the Arnol'd 
diffusion obtained in Jl0| for a classical model. The 
Hamiltonian of this model describes two nonlinear os- 
cillators coupled by the linear term, and governed by an 
external force f(t), 

H=f + ^ + ^ + ^- flX y-xf(t). (1) 

Here p x and p y are momentums in the x and y directions, 
and \i is the coupling constant. The driving term consists 
of two harmonics of the same amplitude /qi 



f(t) = / (cos^i< + cosf2 2 t), 



(2) 



with commensurate frequencies, mf^i = n£l 2 , so that the 
period is T = 2im/Vli — 2nm/Vl 2 . 

Without the coupling (/i = 0) and in the absence of the 
perturbation (/o = 0), the motion of each oscillator is in- 
tegrable and can be found analytically. The quartic form 
of the potentials has been chosen in order to have simple 
analytical expressions in comparison with more realistic 
models with additional quadratic terms x 2 /2 and y 2 /2 in 
the Hamiltonian. 

An interesting feature of the system of quartic oscilla- 
tors is a small contribution of higher harmonics in spite 
of a strong nonlinearity. Indeed, the solution for x(t) has 
the form (for details see ||]), 
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As a result, new actions I\ and I2 are expressed as fol- 
lows, 

h = ^V^, h = (6) 



For the coupling resonance we have I x sa I yi hence, 
I\ -C Iii an d the resonance Hamiltomian H res gets the 
form, 
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with B{I 2 ) = 8/9 • AI 2 Z/A and V(h) = fia 2 (I 2 )/2. 

The average of Eq. (]?]) over the fast phase 9 2 gives the 
"pendulum" Hamiltonian, 



H = 2AI. 
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It defines the frequency Cj of small oscillations of 9\ and 



Co = Vbv = P 



0.85, 



(9) 



2K(1/V2) 

and the half-width Alu of the coupling resonance, 

ALU = (3y%. (10) 

Comparing Eqs.(^) and (^), one can understand that the 
time-dependent perturbation destroys the separatrix of 
the coupling resonance, and gives rise to a chaotic mo- 
tion in the vicinity of the separatrix. Specifically, the 
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action 1% reveals a weak Arnol'd diffusion along the cou- 
pling resonance inside the stochastic layer. 

Thus, the long-term dynamics we are interested in, 
is controlled by three resonances, the coupling and two 
driving ones. The first-order driving resonances are 

determined by the condition lo x (I x ) = Qi, Q2 where 
1/3 

uj x = 4A/3- I x ■ In fact, the coupling resonance u> x = ui y 
serves as a guiding resonance along which the diffusion 
takes place. 

All three resonances are characterized by their widths, 
they can overlap with each other if the coupling constant 
fi and perturbation strength /o are large enough. In or- 
der to observe the Arnol'd diffusion, one needs to avoid 
such an overlap since it leads to a strong global chaos. 
The condition for the overlap of the resonances reads 



where 
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is the half- width of the i-th driving resonance, <5f2 = 
|Oj — £1 2 |- From Eqs.(|To|) and (12) one can obtain for 
the overlap, 
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The Arnol'd diffusion occurs in the case when inequality 
( |l3| ) does not satisfy. 

In principal, the Arnol'd diffusion arises in our model 
even for one driving resonance. However, in this case, 
the rate of the diffusion will be strongly dependent on 
the distance between the position of a trajectory inside 
the stochastic layer, and the driving resonance in the fre- 
quency space. Instead, for two driving resonances the 
Arnol'd diffusion is almost homogeneous if one starts in 
between the two driving resonances. This simplifies the 
analytical treatment of the diffusion, that has been per- 
formed in Refs. j|o|. Leaving aside technical details, we 
briefly comment below the approach used in ]l0| . 

In order to obtain the diffusion coefficient for a diffu- 
sion along the coupling resonance, one needs to find the 
change of the total Hamiltonian over the half-period of 
the unperturbed motion near separatrix. Therefore, the 
diffusion coefficient can be evaluated as follows, 
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where T a is the averaged period of motion within the sep- 
aratrix layer. The change of the Hamiltonian depends on 
an initial phase, however, successive values of the phase 
can be treated as random and independent. The varia- 
tion of the total energy is then determined by the sum 



over many periods for which successive phases can be 
obtained via the separatrix map. Analytical estimates 
obtained in [|o) , give the following expression for the dif- 
fusion coefficient (in action), 



(15) 



Here w s — 47wA 2 e~ 1lA / 2 is the half-width of chaotic layer 
of the coupling resonance, that mainly depends on the 
adiabaticity parameter A = 6fl/2u>, with u determined 
by EqM. 
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FIG. 1. Three sections of the phase space for dif- 

ferent values of I2 are shown for the Hamiltonian (^) with 
the force (|) with /x = 2 ■ 1(T 4 , f = 2 ■ 10" 6 . The Arnol'd 
diffusion takes place along the /2-direction, see details in the 
text. 

We performed numerical study of the classical Arnol'd 
diffusion in the system described by the Hamiltonian (Q). 
The frequencies of the perturbation (||) were chosen as 
ill = 0.2094 and CI2 = 0.2513 resulting in the period 
T = 150. Therefore, u = (fii + fi 2 )/2 = 0.23035 which 
determines the amplitude a sw 0.2719. Correspondingly, 
the initial conditions were taken for the system to be in 
between the two driving resonances. 

To put the system inside the stochastic layer of the 
coupling resonance it is necessary to take 0\ = ±7r. As in 
Ref. jHJ, the ratio fo/n — 0.01 was taken small enough 
in order to avoid the overlapping of three first-order res- 
onances (see Eq.(|l3|)). This also suppresses the influence 
of second-order resonances between unperturbed non- 
linear motion and the external perturbation. 

Schematic structure of the coupling resonance is shown 
in Fig.l. Numerical data are obtained for different ini- 
tial conditions corresponding to the separatrix layer and 
to the resonance region of the coupling resonance. The 
chaotic region inside the resonance is due to second-order 
resonances between non-linear motion of the unperturbed 
Hamiltonian and two-frequency perturbation. The con- 
dition of secondary resonances is n ■ uj' = m ■ SCI. Here £>' 
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is the frequency of oscillations at the coupling resonance 
(near the resonance centre we have £/ — > Co = (i^/Jl) , and 
n, m are integers. One can find that the stochastic region 
inside the resonance corresponds to m = 1 and n — A. As 
one can see, weak diffusion along the coupling resonance 
can occur both within the separatrix layer and inside the 
chaotic region formed by second-order resonances. 
The diffusion coefficient was computed as follows, 
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(AH) 2 
10™ T 



(16) 



Here H is the value of the total Hamiltonian (0), aver- 
aged over time intervals T n of length 10™ -T with n = 2,3. 
The second average in © has been done in the follow- 
ing way. Having the mean value H in each interval T n 
for a fixed n, we computed the difference AH between 
adjacent intervals, and averaged the variance (AH) 2 over 
all differences. This procedure was taken in order to sup- 
press large fluctuations of the energy, and to reveal a 
stochastic character of motion. Specifically, in the case 
of a true diffusion, one should expect D2 ~ D3. 
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FIG. 2. Dependence of the diffusion coefficient on initial 
conditions inside the coupling resonance, see details in the 
text. 

The dependence of diffusion coefficients D2 and D3 on 
initial phase 6\ for I\ = 62 = and 1% — w 3 /(3/? 4 ) at 
t — is shown in Fig.|| The region close to 6\ — — tt cor- 
responds to initial conditions inside the separatrix layer, 
and the interval near 6\ = —tt/2 corresponds to initial 
conditions inside the inner stochastic region (see Fig.0). 
In both these regions Arnol'd diffusion coefficients have 
the same order. Approximate equality D2 ~ D3 indicates 
here that the motion is really diffusion-like |l(| . On the 
other hand, strong difference between D2 and D3 in the 
region \9i\ < ir/4 manifests that the dynamics of the 
system is non-diffusive. 



III. QUANTUM MODEL 

The corresponding quantum model is described by the 
Hamiltonian (compare with Eq.(|l])), 



H = H° + H° - fixy - / a;(cosrM + cos0 2 *)- 
Here 
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(17) 
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and standard relations for momentum and coordinate op- 
erators are assumed, 



[Px, X_ 



= -ifio, 



\Pv, y] = 



(19) 



with the dimensionless Plank constant ho. 

In order to investigate the evolution of the system, 
first, we have to find stationary eigenstates of the un- 
perturbed (/o = 0) system, corresponding to the vicinity 
of the coupling resonance lo x = u> y . At the second stage, 
we will use the Floquet formalism when considering the 
time-periodic perturbation for f ^ 0. Specifically, we 
construct the evolution operator in one period T of the 
perturbation, that allows one to study the dynamics over 
many periods. 



A. Stationary states of the coupling resonance 

It is naturally to represent stationary states of the un- 
perturbed Hamiltonian H s , 



H K = H" 



H° y - fixy 



(20) 



in terms of the eigenstates of uncoupled (ji = 0) nonlinear 
oscillators, 



i>{x, y) = s ^Cn,mi>n(x)i> Q m (y)- 



(21) 



Here tp®(x), ^(j/) are the eigenfunctions of H°, H® 
(which will be calculated numerically), and the coeffi- 
cients c n ^ m satisfy to the following stationary Schrodinger 
equation, 



ECr, 
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n,n' Vm.ra' ^n' ,m' ■> 



(22) 



with E n and E m as eigenvalues of the Hamiltonians H® 
and Hy, respectively. 

At the center of the coupling resonance u no — u mo we 
have n = m Q and h uj no = E' ng , h uj mo = E' mg . Near to 
this resonance it is convenient to expand E n and E m in 
Tailor series up to second order terms. This allows one 
to introduce new indexes p = k + l and k via the relations 
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FIG. 4. Probability |cfe iP | 2 for eigenstates at the coupling 
resonance for q = 0, fj, = 10 -4 and different values of s. a) 
the lowest level (ground state) for s — ; b) near-separatrix 
level, s = 30; c), d) above-separatrix levels, s — ±35. 

Typical structure of eigenstates for different s is shown 
in Fig.[|. Note that ground states in each group cor- 
respond to s = 0, and next stationary states, re- 
ordered according to the energy increase, are labelled by 
1,-1,2,-2,... etc. The eigenstates inside the resonance 
are symmetrical with respect to k = 0. One can see that 
the main maximum corresponds to p = 0, although there 
are small additional maximums at p = ±2, ±4, . . .. The 
degree of derealization in the /c-space for eigenstates in- 
side the coupling resonance strongly depends on the en- 
ergy of eigenstates. Specifically, the closer the energy to 
that corresponding to the separatrix, the more delocal- 
ized is the eigenstate. Above the separatrix all eigen- 
states are characterized by the maximums of probability 
located symmetrically with respect to k = for states s 
and — s, see Fig.||(c),(d). 

B. Evolution matrix 

Now we consider the dynamics of our model in the pres- 
ence of the external two- frequence perturbation acting on 
the x-oscillator. The frequencies f2i and f2 2 are commen- 
surable, so that the external force is periodic with the 
period T = %T X = jT 2 , where T x = 2tt/Oi, T 2 = 2n/Q, 2 
and i, j are integer. The initial conditions were taken for 
the system to be about half-way between the two driving 
resonances, ui = (fli + Sl 2 )/2. 
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FIG. 5. Matrix elements |x 9lS;9 ±i lS ' | that define the tran- 
sition probability along the coupling resonance (for the same 
parameters as in Fig.j3|). Different values of the matrix el- 
ements are shown by grayscaled symbols. Black squares at 
the diagonal of the matrix have the values six orders of mag- 
nitude higher than those of matrix elements labeled by grey 
rhombuses. 

The matrix elements £g,s;g±i,s' m Eq.(pl|) correspond 
to transitions between the states s and s' from neighbor 
groups q and q ± 1. Figure || illustrates relative ampli- 
tudes of the matrix elements |a?o,s ; i.s' I ■ hi accordance 
with our numeration of the states, matrix elements at 
the center of Fig.^ correspond to transitions between the 
lowest states in each group, and matrix elements at the 
corners define the transitions between the states above 
accumulation points. The latter elements quickly de- 
crease with an increase of the difference \s — s'\. 

The "cross" at the center of Fig.[| where matrix ele- 
ments are relatively large, corresponds to a transition be- 
tween separatrix states. The important point is that the 
transition between such states of neighbor groups (along 
the coupling resonance) is much stronger than those be- 
tween other states. This phenomenon is analogous to 
the quantum diffusion inside a separatrix, which was ob- 
served in jl9| for a degenerate Hamiltonian system. 

As a result, numerical procedure for computing the 
dynamics of our model is as follows. First, we solve 



G 



Eqs.([3l|) and construct the evolution matrix U q . s - q '. s '(T) 
by making use of Eq.(|3(]). Then, direct diagonalization 
of this matrix yields the eigenvalues eq and the eigen- 
vectors A® s , Once the eigenvalues eq and eigenvectors 
A^ s are obtained, one gets the evolution operator for one 
period, 



U, 



q,s;q ,s 



exp 



h 



(32) 



By raising U q ^ s - q i , S '(T) to the iV-th power and using the 
ortogonality condition for the eigenvectors A® s , one can 
obtain the evolution operator that propagates the system 
over TV periods of the external perturbation, 



U, 



<3 



ie Q NT 



(33) 



C. Numerical data 

As was shown above, in our model the Arnol'd diffu- 
sion occurs along the coupling resonance, or, the same, in 
q-space. This means that a wave packet initially localized 
at q = 0, spreads diffusively in time. In order to observe 
this dynamics, below we introduce specific variables that 
characterize global structure of wave packets. But, first, 
we discuss the structure of QE eigenstates (in terms of 
these new variables) since it helps us to understand the 
mechanism of quantum Arnol'd diffusion. 

In the g-space each of QE functions A® s can be glob- 
ally characterized by the "mean position" q and by the 
variance o~ q determined as follows p7|JT8[ , 

? = E ?E i^i 2 - < = £(« - ^ E \ A t\ 2 - 

q s q s 

(34) 

Then, it is convenient to plot q versus a q for all eigen- 
functions, see Fig.[| For small values of the coupling /z, 
see Fig.^|(a), the QE functions have a very small variance 
that means a strong localization in the g-space. 

More details are seen in Fig.^(b) which shows a magni- 
fied fragment of Fig.||(a). One can see different groups of 
QE functions, characterizing specific relations between q 
and a q . First group consists of those eigenstates whose 
energies are close to the ground state (with a very small 
variance, a q < ICR 3 ). Another group is characterized 
by an irregular dependence of q on a q (region (1) in 
Fig.^(b)). The states which belong to this group are 
chaotic separtrix eigenstates. Regular dependences q{cr q ) 
for <Jq > ICR 3 correspond to under-separatrix states with 
s > and s < 0. For very large s this regular structure is 
destroyed due to the influence of two driving resonances 
(not shown in Fig.|^(b)). Irregular spread of points in 



the region (2) reflects chaos in the inner region of the 
coupling resonance, that arises due to two secondary res- 
onances hodfl/2 = 3E" (s + 1/2). With an increase of 
H (see Fig.||(c)), the regular structure of QE functions 
disappears. This means that many of eigenstates are 
chaotic. However, the variance at remains limited, thus, 
demonstrating that in the q direction the eigenstates are 
localized. The fact that many points in Fig.^(c) are dis- 
tributed incidentally, should be treated as the manifesta- 
tion of quantum chaos. 
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FIG. 6. Relation between q and a q for the QE functions 
A® 3 in the region q — 0, ±1 for different coupling constants: 
a,b) n = 3 • 10 -5 , c) fj, — 10 -4 . Each point corresponds to a 
specific QE function. FigM b) shows the scaled-up fragment 
of the central part of FigMa). 

Now we discuss numerical results for the dynamics of 
our model. The evolution of any initial state C q '. s ' (0) can 
be computed using the evolution matrix U qiS . q / tS /(NT), 



C q . s {NT) = U q , s;q ,, s >{NT)C q ,, s ,(0). 



(35) 



q',s' 



We repeat again that our numerical data refer to the 
regime when the values of /i and /o are small enough, so 
that main three resonances are not overlapped. 

Quantum dynamics for different initial conditions is 
shown in Fig.j^. Here we show typical dependencies of 
the variance of the energy (AH) 2 = Ti 2 ) uj 2 /S. q in nor- 
malized units, versus time t measured in the number 
N of periods of the external perturbation. The quan- 
tity A q is defined similar to that for the QE eigenstates, 
A 9 =£,(«" q? Es \C q ,s\' 2 where q = £ 9 g£ s |C 9 , S | 2 . 

The data clearly show a different character of the evo- 
lution for three initial states taken from below and above 
the separatrix, as well as from the separatrix layer. For 
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the states taken from the center of the resonance and 
well above the separatrix, the variance (AH) 2 quasi- 
periodically oscillates, in contrast with the separatrix 
state. In the latter case, after a short stochastization time 
the variance of the energy increases linearly in time, thus 
manifesting a diffusion-like spread of the wave packet. 



FIG. 8. Time dependence of the variance of the energy for 
separatrix initial states and different values of coupling pa- 




1000 



FIG. 7. Normalized variance A 9 of the energy versus the 
rescaled time N for 3 different initial states (for fi = 1.25- 10 -4 
and fo = 1.25T0" 6 ). Curves 1-3 correspond to an initial state 
near the center of the coupling resonance, above the separa- 
trix, and from the separatrix layer, respectively. 

More results are presented in Fig.|| where a diffusion- 
like increase of the energy is shown for separatrix initial 
states and different values of /i. One can see that linear 
increase of the variance A q (N) is typical, and it occurs af- 
ter a short time which is associated with the time of a fast 
spread of packet over the separatrix layer in transverse 
direction. Such a behaviour is typical for the classical 
Arnol'd diffusion. The data allows one to determine the 
diffusion coefficient as D = (AH) 2 (N)/(N • T) , by mak- 
ing use the fit to a linear dependence A q (N) for N > 50. 
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We have calculated separately quantum and classical 
diffusion coefficients and found that the quantum Arnol'd 
diffusion roughly corresponds to the classical one. How- 
ever, the data clearly indicate that the quantum diffusion 
is systematically weaker than the classical Arnol'd diffu- 
sion, see FigJ|. 



Log D 




75 100 

FIG. 9. Quantum (squares) and classical (solid line) diffu- 
sion coefficients for different values of \i. 

One should stress that the quantum Arnol'd diffusion 
takes places only in the case when the number M s of en- 
ergy stationary states in the separatrix layer is relatively 
large. For the first time, this point was noted by Shuryak 
Ip^l who studied the quantum-classical correspondence 
for nonlinear resonances. In this connection we have to 
estimate the number of the energy states that occupy the 
separatrix layer. 

However, it is a problem to make reliable analytical es- 
timates of the width of the separatrix layer for our model. 
As was shown in pcj , in the case of many-frequency per- 
turbation, secondary resonances play the dominant role 
in the formation of chaos in a separatrix layer. Indeed, 
as was shown above, the role of such resonances is quite 
strong in our case of the two-frequency perturbation. For 
this reason, we have performed a direct numerical calcu- 
lation of the width of stochastic separatrix layer for the 
classical Hamiltonian, and used these results in order to 
estimate the number of separatrix energy levels in the 
corresponding energy interval. 

We have found that for fi > 1.25 • 10~ 4 the number M s 
of stationary states in the separatrix chaotic layer is more 
than 10, therefore, one can speak about a kind of stochas- 
tization in this region. On the other hand, with a de- 
crease of the coupling (see data in Fig.^for 1/y/ji > 100), 
the number M s decreases and for fx « 3 ■ 10~ 5 it is of the 
order one. For this reason the last right point in FigJ9 cor- 
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responds to the situation when classical chaotic motion 
along the coupling resonance is completely suppressed by 
quantum effects (the so-called "Shuryak border" JT^j). 



FIG. 10. Dynamical localization of wave packets inside the 
separatrix layer of the coupling resonance. The normalized 
variance of the energy is shown for different values of the cou- 
pling parameter: (1) fj, = 1.5 • 1(T 4 , (2) fi = 2.25 ■ 1(T 4 . 
Horizontal dashed lines indicate the mean value A 9 . 



D. Dynamical localization 



Since the diffusive motion along the coupling resonance 
is effectively one-dimensional, one can naturally expect 
the Anderson-like localization. We have already noted 
that the variance of QE eigenstates of the evolution oper- 
ator is finite in the g-space. This means that eigenstates 
are localized, and the wave packet dynamics in this direc- 
tion has to reveal the saturation of the diffusion. More 
specifically, we expect that the linear increase of the vari- 
ance of the energy ceases, after some characteristic time. 

This effect known as the dynamical localization, has 
been discovered in |2l],|2^] for the kicked rotor, and was 
studied later in different physical models (see, for exam- 
ple, |l| and references therein). One should note that the 
dynamical localization is, in principle, different from the 
Anderson localization, since the latter occurs for mod- 
els with random potentials. In contrast, the dynamical 
localization happens in dynamical (without any random- 
ness) systems, and is due to interplay between (week) 
classical diffusion and (strong) quantum effects. 

In order to observe the dynamical localization in our 
model (along the coupling resonance inside the separatrix 
layer), one needs to study long-time dynamics of wave 
packets. Our numerical study for large times N « 10 4 
have revealed that after some time t ~ to as 10 3 • T, 
the diffusionlike evolution stops for all range of coupling 
parameter fi. Instead, for larger times, the variance A q 
starts to oscillate around the mean value A q . Period 
of such oscillations depends on \i non-monotonically and 
vary from 10 3 • T to 10 4 • T. Two examples of such a 
long-time dynamics are given in Fig.fLOl 



One can argue that the localization length is of the or- 
der of the width of wave packet after the saturation of 
the diffusion |p2| , [22|| . Therefore, the localization length 
l s can be associated with the square root of A q . Wc 
have numerically found the exponential dependence A q 
on 1/s/Jl, see Fig.pl]. This result is not surprising be- 
cause log D is proportional to 1/ ^JJi (see Fig.||) and to is 
practically independent on fi. 



10* 
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FIG. 11. The dependence of A q on 1/^/Jl. 

It should be noted that the spread of packets in k and I 
space occurs even without the coupling between two res- 
onances, due to the influence of the time-dependent per- 
turbation (see Section III- A). For this reason one should 
compare this spread to that determined by the Arnol'd 
diffusion. Our additional study clearly show that these 
two effects are very different, see Fig.|l2[ 
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FIG. 12. Distribution of the probability W q , averaged over 
time t — 10 4 ■ T, versus the group number q: (1) /i = 0, 
fo = 10" 6 ; (2) p = 0, f = 2 • lCT 6 ; (3) fi = lO" 4 , /„ = 10" 6 ; 
(4) M = 2-10- 4 , /o = 2-10" 6 . 

Here, we plotted the profile W q = J2 S \Cq,s\ 2 of the 
wave packet after the saturation, versus the group num- 
ber q. This hgure demonstrates the main effect of the 
Arnol'd diffusion along the coupling resonance. One can 
see that in all cases there is an exponential localiza- 
tion of packets in the g-space. This allows to introduce 
the localization length defined from the decrease of the 
probability in the tails of packets. The data illustrate a 
strong increase of the localization length in the presence 
of the coupling between two oscillators, in comparison 
with the case of completely independent oscillators for 
fi = (q = k + I and I = 0). 

IV. SUMMARY 

We have studied the Hamiltonian system of two cou- 
pled nonlinear oscillators, one of which is under the influ- 
ence of the time-dependent perturbation with two com- 
mensurate frequencies. In the classical description, the 
separatrix of the coupling resonance is destroyed due to 
the perturbation, and the Arnol'd diffusion occurs along 
this resonance inside a narrow stochastic layer. Our nu- 
merical data performed for the quantum analog of the 
system, allow us to make the following conclusions about 
the properties of the quantum Arnol'd diffusion. 

By studying the quasienergy eigenstates of the evo- 
lution operator, we have found an irregular structure 
of those eigenstates which correspond to the stochastic 
layer. These eigenstates turned out to be exponentially 
localized along the coupling resonance, with the localiza- 
tion length, strongly enhanced in comparison with the 
case of non-coupled oscillators. 

The study of wave packet dynamics for different ini- 
tial states have revealed a diffusion-like spread of packets 
along the coupling resonance, if initial states correspond 
to the stochastic layer. We have found that the depen- 
dence of the diffusion coefficient on model parameters, 
roughly follows the classical dependence. However, the 
quantum diffusion is systematically slower than the clas- 
sical one. This fact is obviously due to the influence of 
quantum effects. 

It should be stressed that the quantum Arnol'd diffu- 
sion occurs in a deep semiclassical region, specifically, for 
the case when the number M s of chaotic eigenstates in- 
side the stohastic layer is sufficiently large (of the order 
of 10 or larger). With a decrease of the coupling param- 
eter, the diffusion coefficient strongly decreases, and for 
M s < 1 the diffusion disappears. Therefore, we can see 
how strong quantum effects destroy the diffusive dynam- 
ics of wave packets. 



Another manifestation of quantum effects is the dy- 
namical localization that persists even for large M s . 
Specifically, we have observed that the quantum diffu- 
sion occurs only for finite (although large) times. On 
a larger time scale the diffusion ceases, and after some 
characteristic time it terminates. This effect is similar to 
that discovered in the kicked rotor model plj and stud- 
ied later in other physical systems (see, for example, jlj] 
and references therein). However, in our case the dy- 
namical localization arises for a weak chaos inside the 
separatrix layer, in contrast to previous models with a 
strong (global) chaos in the classical description. 

Our results may find a confirmation in experiments on 
one-electron dynamics in 2D semiconductor quantum bil- 
liards, where the charged particle motion is determined 
by the Hamiltonian of the type (|l7|), (|T§|). It is also possi- 
ble that the quantum Arnol'd diffusion occurs in nuclear 
dynamics of complex molecules, driven by laser fields Q . 
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